Massless Fermions in multilayer graphitic systems with misoriented layers. 
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We examine liow tlie misorientation of a few stacked graphene layers affects the electronic structure 
of carbon nanosystems. We present ah initio calculations on bi- and trilayer systems to demonstrate 
that the massless Fermion behavior typical of single layered graphene is also found in incommen- 
surate multilayered graphitic systems. We also investigate the consequences of this property on 
experimental fingerprints, such as Raman spectroscopy and scanning tunneling microscopy (STM). 
Our simulations reveal that STM images of turbostratic few layer graphite are sensitive to the layer 
arrangement. We also predict that resonant raman signal of graphitic samples are more sensitive to 
the orientation of the layers than to their number. 
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The electronic properties of 2D graphene and 3D 
graphite have been extensively studied for more than 50 
years It is fascinating that the extraordinary properties 
of carbon nanotubes^. have been deduced from those of 
2D graphene many years before macroscopic samples of 
very thin Few Layer Graphite (FLG) could be obtained in 
the laboratory'^. A particular interest has been recently 
given to Single Layer Graphene (SLG) because of their 
massless Fermion behavior, the Vb dependence of the 
Landau levels'''^, and the observation of abnormal Quan- 
tum Hall Effect (QHE), even at room temperature^. 

In that context, a precise investigation of the layer to 
layer interaction on the existence of massless Fermion 
carriers is of paramount importance. For 3D graphite, 
the most stable Bernal phase (AB stacking) as well as the 
rhombohedral (ABC stacking) have been proven to show 
complex electron and hole bands near the Fermi level 
rather than linear, massless fermion ones, due to inter- 
layer interaction^. We have recently shown that regular 
(AB or ABC) stackings also break the linear character 
of the dispersion of electronic bands for FLG with 2 to 
4 layers and that ambipolar electronic conduction could 
be related to AB stacked FLG^. 

In this Letter, we present an electronic structure analy- 
sis of misoriented (turbostratic) 2 and 3 layer FLG. A re- 
cent surface X-rays analysis of multilayer graphene grown 
on SiC shows that misorientation is plausible^. We show 
here that the linear dispersion of SLG is preserved in tur- 
bostratic multilayer systems despite the presence of adja- 
cent layers. It follows that massless Fermion carriers are 
predicted for disoriented multilayer systems. These find- 
ings raise the question of the interpretation of the experi- 
mental observations performed on FLG samples and chal- 
lenge the direct relation between SLGs and Dirac mass- 
less Fermions. More generally, the electronic properties 
(and consequently the optical, vibrational, and transport 
properties) of a given FLG film is found to be controlled 
mainly by the (mis)oricntation of the successive layers 
rather than the number of layers. In the present work , 
we discuss also the implications of this possible misorien- 
tation on experimental signatures, in particular on STM 



and Raman fingerprints. 

The main technical difficulty when modeling tur- 
bostratic structures, is to combine their incommensurate 
character and the necessary finiteness of the supercell in 
solid state calculations. An elegant solution, proposed by 
Kolmogorov and Crespi^- consists in the definition of a 
graphene hexagonal supercell without mirror symmetry 
(except the basal plane). The smallest supercell of this 
type contains 14 atoms, and is defined with its basis vec- 
tors Ai = 2ai +aL2 (denoted as (2, 1) supercell hereafter) 
and A2 = — ai -I- 3a2, as shown on Fig[TJ The super- 
cell (3, 1) possesses exactly the same basis vectors, but 
is rotated by the angle cos"^(ll/14) ~ 38.21°. When 
stacking these two supercells, we obtain a bilayer struc- 
ture with short-range incommensurability (FigU]). The 
(7, 0) and (5, 3) supercells also present exactly the same 
size but a larger number of atoms (98 C atoms per layer), 
and can be assembled onto bi- or trilayer compounds ful- 
filling the translational symmetry requirements^^. 

In this Letter, we focus our attention on bi- and tri- 
layer FLG structures based on the (2, 1) and (7, 0)/(5, 3) 
supercells. Our aim is to highlight significant differences 
in the band structures of commensurate and incommen- 
surate graphitic systems. All the calculations were per- 
formed with density functional theory (DFT), within the 
local density approximation (LDA) schemei^, and using 





turbostratic AA' supercell Bernal AB supercell 



FIG. 1: a) The non-symmetric supercell (2,1) (light grey), 
and its two primitive vectors Ai and A2. The (3, 1) cell is 
drawn in dark grey, b) A quasi incommensurate bilayer AA' 
structure, made by stacking the (2, 1) on the (3, 1) supercell. 
c) An AB bilayer structure supercell. 
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turbostratic AA' superccll 




FIG. 2: The 2D electronic band structures of Bernal AB (a) 
and turbostratic AA' (b) bilayer graphene. Both are based on 
the supercell represented in Fig. 1. The hexagonal Brillouin 
zone of the system is shown in the inset. 



norm-conserving pseudopotentials^^. According to the 
number of atoms involved in the supercell, the eigen- 
problem was expanded on plane waves (PWs, using the 
code ABINIT^'*) or localized pseudoatomic orbitals (PAOs, 
using the code SIESTA^^). In the first case, plane waves 
cutoff energy was set to 35 Hartree. In the second case, 
the basis was composed of atom-centered double-C func- 
tions. 

We carried out a structural optimization of each bilayer 
(2, l)-based supercell (PW calculations with tolerance for 
the forces set to 10~^ Hartree/Bohr and a 12x12x1 grid 
sampling of the Brillouin zone), prior to the analysis of 
their electronic structure. No significant atomic rear- 
rangement is observed within a plane, the only structural 
reorganization taking place in the average inter-layer dis- 
tance (3.33 A for the AB, and 3.42 A for the turbostratic 
stacked bilayer). The corresponding band structures of 
the two (2,l)-based bilayers are plotted on Fig. O We 
have carefully verified that PAO basis yields identical re- 
sults. 

Fig. [5^ reproduces the well-known results for AB bi- 
layeriSii^ within a supercell \/l x \/l larger then the primi- 
tive cell and is displayed for comparison. As in previously 
published calculations, we notice a ^ 0.8 eV splitting of 
the electronic bands, related to the layer-layer interac- 
tions in the AB geometry. The bilayer graphene with a 
38.21° misorientation (called AA' hereafter) presents a 
totally different signature near the Fermi level (Fig[2lD). 
Surprisingly, the band structure here is found to be simi- 
lar to the one of SLG but doubly degenerate. The signif- 
icance of this result is striking, since a massless Fermion 
character is found for misoriented bilayer systems with 
a Fermi velocity of ^^9.6x10^ m s~^, strictly identical to 
the value deduced for SLG within the same formalism, 
and similar to experimental values§-i2?. We have checked 



that this unexpected behavior is also found in the tur- 
bostratic (5, 3)on(8, 3) bilayer that presents a misorienta- 
tion angle of 43.57°. The apparent absence of interaction 
between layers at the Fermi level can be related to the 
loss of short range corrugation between successive layers. 

Turning to the trilayer FLG system, we use the inter- 
layer distance found for bilayer cases, depending of the 
stacking geometry (AB or turbostratic). The PAO elec- 
tronic calculations shown here have been checked against 
PW calculations. On FiglU we present the results for a 
purely AA'A" turbostratic (5, 3)on(8, 3)on(7, 0) trilayer 
with two respective misorientations of 43.57° and 38.21°, 
and for a mixed case (5, 3)on(7, 0)on(7, 0), made from a 
bilayer with AB stacking and a third one misoriented 
(38.21°), called ABA'. The results for the calculation for 
ABA or ABC stacking can be found in references^^-'-^^. 

The trilayer with a mixed stacking ABA' displays elec- 
tronic bands (Fig. [31d) similar to a superposition of a AB 
bilayer (Fig. [2^) and a single-layer graphene (linear dis- 
persion). The band crossing of the graphene- like bands 
lies 12 meV above the Fermi level, indicating a small 
charge transfer of ~ 10^^ e cm~^. The other (quadratic) 
bands present a gap opening that is consistent with the 
behavior found for AB bilayer graphene under an elec- 
tric field^^. We have checked that the same systems 
but with the A' layer away from its equilibrium position 
presents a vanishing small charge transfer, i.e. a pro- 
gressive closure of the 'gap' of the quadratic bands and a 
Dirac point of the linear bands moving towards the Fermi 
level. Interestingly, the turbostratic AA'A" trilayer fol- 
lows the trends observed for the bilayer and displays lin- 
early dispersive bands and massless Fermion behavior. 
The (doubly degenerate) bands associated with the ex- 
ternal layers are slightly shifted upward (~12 meV) while 
the bands related to the central layer are slightly shifted 
downward(~25 meV) . This is also the signature of a small 
electron transfer, from the external layers to the central 
one, consistent with the previous ABA' case. We have 
also checked that the deviation from the linear dispersion 
of the TT-bands for \E\ > 0.5 eV is mainly due to the trig- 
onal warping effect of graphene and not to an interlayer 
effect. 

The results presented above demonstrate that tur- 
bostratic stacking of graphene sheets leads to a linear 
dispersion of the electronic levels similar to the behavior 
of SLG. Furthermore, the mixed stacking case shows a 
superposition of the two types of carriers, i.e. a massless 
Fermion and a massive, normal electronic behavior. Con- 
sequently, direct experimental evidence of Dirac Fermion 
behavior can no longer be considered as a discriminating 
property between single and multilayered systems. In 
particular, the number of tt bands in Angle resolved pho- 
toemission spectroscopy (ARPES) measurement^ is not 
a definitive evidence of the number of layers but should 
be cross-checked with other experimental techniques. 

The absence of a net effect of interlayer interaction on 
the electronic dispersion of turbostratic FLG also sug- 
gests that the ^/B dependence of the Landau levels will 
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also occur for multilayer structures and that the abnor- 
mal QHE cannot be excluded. Remarkably, this conse- 
quence of our study is supported by recent experimental 
data. First, ^/~B Landau level separation have been ob- 
served in 3-5 graphene layers^. Second, Dirac Fermion 
behavior, integer QHE-^ and infrared probe anomalous 
magnetotransport have been also reported in a Highly 
Ordered Pyrolytic Graphite (HOPG) bulk 3D sample^°. 
Third, the unconventional QHE has been reported on bi- 
layer grapheneSi, challenging the common interpretation 
in terms of AB stacking. Finally, other studies^ relate 
the coexistence of two types of carriers (Dirac Fermion 
and normal electron) in HOPG bulk graphite, consistent 
with the findings presented herein. 

Recently, the observation of a single 2D band (~ 2700 
cm~^) in Raman spectroscopy has been proposed as an 
experimental tool for the determination of the number of 
layers in FLG samplesi^i^i^. For two-layer (and thicker) 
films, the interlayer interaction induces a splitting of the 
electronic bands. This means that in a double resonance 
Raman process, the resulting effect is a splitting of this 
2D peaki^i^i. Our results show that the effect of the 
interlayer interaction on the splitting of the electronic 
bands crucially depends on the relative orientation of 
the layers. For instance, since a turbostratic bilayer pos- 
sesses degenerate electronic bands near the Fermi level, 
the 2D mode of this structure gives rise to a single peak 
that is indistinguishable from the signature of SLG sys- 
tems. Similarly, a trilayer turbostratic structure will also 




FIG. 3: Representation of the supercells and the correpond- 
ing electronic band structure for trilayer graphites. The tur- 
bostratic AA'A" system (43.57° and 38.21° misorientation 
angles) is shown in (a), and the mixed ABA' structure is 
shown in (b). 



present a single 2D band. Indeed, as we have discussed 
above, the splitting of the electronic bands in that case 
is due to a charge transfer and no evidence of covalent 
mixing has been found. The 2D peak of the mixed ABA' 
trilayer will present three main features. Two of them are 
associated with the AB-like bands, the third one with the 
A' layer, identical to the SLG case. 

Interestingly, the splitting of the 2D band of graphitic 
materials has been found to be a good criteria of crys- 
tallinity perpendicular to the basal planes. For instance, 
in Ref. [23, the authors have correlated the interlayer 
distance in bulk graphite with the doubling of the 2D 
band and they found that for d > 0.338 nm a single 2D 
band is observed where a double 2D band is observed for 
d < 0.338 nm. From these findings and the results pre- 
sented here, we can propose that HOPG graphite with 
disorder in the staking (and large d) will present a single 
2D line whereas Bernal graphite with (partial) AB order 
will present a doubling of the 2D line. ^From this dis- 
cussion, it is clear that a unequivocal determination of 
the number of layers is not possible from Raman spec- 
troscopy alone. 

We now analyse the STM images of the trilayer sys- 
tems studied above. Recent experimental data on mono- 
and bilayer graphite on SiC substrate show a a/3 x \/3 
modulation of the STM intensity. This effect has been as- 
sociated with the substrate reconstruction^LSliS^ rather 
than to the relative arrangement of the graphitic layers. 
Here, we have computed the STM images of trilayer sys- 
tems without the substrate in order to discriminate the 
intrinsic response of FLG from the effect of interaction 
with the substrate. This approach is made possible by 
the prerequisite of Tersoff-Hamann theory since the im- 
age formation is governerd by integration of local density 
of states, and does not require a closed circuit to simulate 
the curren1>2S. Here we made use of the BSKAN package 
implementation of Tersoff-Hamann^Ml. Constant cur- 
rent images for systems represented in Fig. [3] are shown 
in Fig. m for positive and negative tip polarities. For sys- 
tem ABA' the two external faces are not equivalent and 
the corresponding images are given for both faces. At 
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ABA' (A' face) ABA' (A face) 



FIG. 4: (color online) ^From left to right: constant cur- 
rent STM image simulation of turbostratic AA'A" and mixed 
ABA' systems. (Un)occupied states are images at (positive) 
negative tip potential. The color scale spans the corrugation 
of each system: 1.6 (1.4) A, 2.3 (2.4) A, and 1.8 (2.1) A for 
each system at positive (negative) bias. 
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low bias (100 meV was used here), the images reflect the 
property of the states close to the Fermi level. For sys- 
tem AA'A" ', both positive and negative bias images look 
the same, with all the atoms being imaged, just as in 
SLG. The situation is quite different for ABA' systems. 
The image for occupied states of the A' face shows bright 
spots on every other atom, while a complete honeycomb 
lattice is imaged for unoccupied states. The situation 
is even more complicated for images computed on the 
A face of the ABA' system. While every other atom 
is now revealed at tip bias corresponding to unoccupied 
states, the occupied state image clearly shows a small- 
scale Moire-like pattern. It is important to note that the 
images have been simulated for zero temperature con- 
ditions: given the very small energy difference between 
the states responsibles for the images at low bias, the 
room temperature image will most probably look like an 
average of the images reproduced here. 

In summary, we have performed ab initio calculations 



of the electronic band structures of misoriented (tur- 
bostratic) FLG. We surprisingly found that, as a con- 
sequence of the incommensurability of the layers, such 
multilayer systems possess a massless Fermion carriers 
property. Since that property was up to now consid- 
ered as a unique signature of single layer graphene, the 
present findings challenge the interpretation of the exper- 
imental data on FLG. They also release the restriction to 
SLG sample to build devices based on massless fermion 
physics. 
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Namur. S.L. and L.H. acknowledge the financial support 
from the Belgian FNRS. V.M. acknowledges the Division 
of Materials Sciences, US Department of Energy, under 
Contract No. DEAAC05-G0OR22725 with UT Battelle, 
LLC at ORNL. 
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